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QPOs during magnetar flares are not driven by mechanical 
normal modes of the crust. 
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ABSTRACT 

Quasi-Periodic Oscillations (QPOs) have been observed during three powerful magne- 
tar flares, from SGR0526-66, SGR1806-20 and SGR1900+14. These QPOs have been 
commonly interpreted as being driven by the mechanical modes of the magnetar's 
solid crust which are excited during the flare. Here we show that this interpretation is 
in sharp contradiction with the conventional magnetar model. Firstly, we show that a 
magnetar crustal mode decays on the timescale of at most a second due to the emis- 
sion of Alfven waves into the neutron-star interior. A possible modification is then to 
assume that the QPOs are associated with the magnetars' global modes. However, we 
argue that at the frequencies of the observed QPOs, the neutron-star core is likely to 
support a continuum of Magneto-Hydrodynamic (MHD) normal modes. We demon- 
strate this on a completely solvable toy model which captures the essential physics of 
the system. We then show that the frequency of the global mode of the whole star is 
likely to have a significant imaginary component, and its amplitude is likely to decay 
on a short timescale. This is not observed. Thus we conclude that either (i) the origin 
of the QPO is in the magnetar's magnetosphere, or (ii) the magnetic field has a special 
configuration: either it is expelled from the magentar's core prior to the flares, or it's 
poloidal component has very small coherence length. 



1 INTRODUCTION 



In a prophetic paper, Duncan and Thompson (1992) argued 
that magnetars — a class of neutron stars with super-strong 
(10 14 — 10 15 G) magnetic fields — must exist. Treated at first 
with some skepticism, the magnetar paradigm has proved 
extremely successful in explaining the rich phenomenology 
of Anomalous X-ray Pulsars (AXPs) and Soft Gamma-ray 
Repeaters (SGRs) (e.g., Thompson and Duncan 1995, TD). 
In particular one phenomenon which has an appealing ex- 
planation within the magnetar paradigm are Giant Flares 
from SGRs. They are thought to be powered by an im- 
pulsive release of magnetic energy stored in the neutron 
star. The release may be triggered by a fracture in the 
magnetically-stressed crust (TD) or by a sudden reconnec- 
tion in a twisted magnetosphere (Lyutikov 2003). So far 3 
Giant Flares have been observed: from SGR0526-66 (Mazets 
et al. 1979), SGR1900+14 (Hurley et al 1999, Feroci et 
al. 1999), and SGR1806-20 (Hurley et al. 2005, Palmer et 
al. 2005). In each of the flares' light-curves, there is intrigu- 
ing evidence for QPOs with frequencies of tens of Hertz. 
Firstly, the 43.5Hz QPO was reported by Barat et al. (1983) 
in the 1979 flare from SGR0526-66. Then recently Israel 
et al. (2005) have found several QPOs in the 2004 Giant 
Flare from SGR1806-20, with frequencies from 18 to 90 Hz. 
These QPOs were detected with very high signal-to-noise 
and lasted for ~ 100 seconds; their presence was recently 
confirmed by Watts and Strohmayer (2005, WT) who used 



data from a different satellite. WT have also detected a rel- 
atively high-frequency QPO, of 625Hz. Lastly, prompted by 
the Israel et al. observations, Strohmayer and Watts (2005) 
have found multiple QPO in the light-curve of the 1998 Gi- 
ant Flare from SGR1900+14. 

The QPOs have been widely associated with torsional 
modes of the neutron-star crust [Barat et al. (1983) made 
the first suggestion, followed by a more detailed analysis of 
Duncan (1998), Israel et al. (2005), Strohmayer and Watts 
(2005) and Piro (2005)]. Indeed, it is attractive to associate 
a stable QPO with a mechanical mode of a neutron star, and 
the frequencies of crustal torsional modes are of the right or- 
der of magnitude* However, we argue below that this inter- 
pretation is not viable within the magnetar paradigm. In the 
next section, we show that the crustal mode of a magnetar 
is not stable but decays within seconds by emitting Alfven 
waves into the core. This in direct contradiction with the ob- 
servations. Chris Thompson (private communications) has 
suggested that one may associate the QPOs with the MHD- 



The 625Hz QPO found by WT is already somewhat problem- 
atic for the crustal-mode picture. One can associate it with a 
crustal shear mode which has 1 radial node (n=l; see e.g. Piro 
2005). However, a large multitude of / < 5, n = 1 shear modes 
exist around that frequency, with spacings of several Hz. There is 
no physical reason why just one of them should be excited, and 
not many. 
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elastic modes of the whole neutron star. However, in section 
3 we argue that the frequency of such global mode is likely 
to have an imaginary component, and the mode decays on 
a short timescale. This is due to the continuous spectrum of 
MHD modes supported by the liquid core. 



2 DECAY OF MAGNETAR CRUSTAL MODES. 

Magnetic fields mechanically connect the elastic crust with 
the liquid core. Both media can be considered as perfect 
conductors, and ideal MHD is applicable with magnetic field 
lines frozen into the media. The Alfven wave speed in the 
core is given by 



-BBcr 
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3 x 10' 



cm/s. 



(1) 



Here B = 10 B15G is the magnetic field inside the core, 
p = 10 14 pi4g/cm 3 is the density, B crit ~ 10 15 G is the field 
inside superconducting fluxtubes, and r? is the fraction of the 
neutron-star core involved in the Alfven-wave motion. If the 
neutrons are superfluid, then they decouple from the Alfven 
wave and r\ ~ 0.1. If protons do not form a superconducting 
fluid, -B cr it should be substituted by B in the above equation. 

The timescale for an Alfven wave to cross the neutron- 
star core is ~ 0.05s. Therefore, even if the oscillation was 
initially localized in the crust, the rest of the core will get 
involved in less than a second. It is instructive to estimate 
the power W radiated in the Alfven waves from the crust: 



W ~ pu\ 2 v a A, 



(2) 



where co is the angular frequency of the oscillation, £ is 
the amplitude of the crustal oscillation at the crust-core in- 
terface, and A is the inner surface area of the crust. The 
timescale on which the energy is drained from the crustal 
mode is given by 

M 



0.01s. 



(3) 



pV a A 

Here M is the mass of the crust. 

In a completely rigorous treatment of the problem var- 
ious corrections of order 1 would occur. Partial reflection 
of the Alfven waves is possible, but it will not be a ma- 
jor correction since the wavelength and the relevant scale- 
heights (density, magnetic pressure) are comparable to the 
neutron-star radius, and thus there is no impedance mis- 
match which is needed for significant reflection. We note 
that the situation here is qualitatively different from that 
considered by Blaes et al. (1989), who considered small- 
wavelength shear waves propagating from the depth of the 
crust into the neutron-star magnetosphere, and found large 
reflection coefficients. Here we deal with global modes of 
the crust, which feature almost radius-independent horizon- 
tal displacements. The motion of the inner surface of the 
crust, together with frozen-in magnetic field-lines, should be 
considered as boundary condition for launching the Alfven 
waves into the core. 

The geometry of the magnetic field is another obvious 
factor, but for a generic configuration we do not expect a 
qualitative change in radiated Alfven-wave power. If how- 
ever magnetic field was confined to the crust or had a very 
small coherence length in its poloidal component, then the 
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Figure 1. Toy model of the magnetar core. Perfectly conducting 
incompressible fluid is sandwiched between perfectly conducting 
top and bottom plates. Magnetic field — shown by thin arrows — 
is vertically directed and is a function of y only. The z-axis is 
perpendicular to the page. 



radiated Alfven-wave power could decrease significantly. We 
do not know any compelling argument for either. Thus we 
conclude that even if the crustal mode was originally excited, 
magnetic stresses would significantly reduce the mode am- 
plitude in a fraction of a second and redistribute the energy 
within a liquid core. Therefore, barring special magnetic- 
field geometry, a pure crustal mode cannot be associated 
with a stable QPO during the flare. One may hope to as- 
sociate the global MHD-elastic mode with the QPO. This 
proposal (due to Chris Thompson) runs into potential diffi- 
culties which we discuss in the next section. 



3 GLOBAL MODES. 

Computing global MHD modes of a star is clearly a difficult 
task. In recent tour-de-force study Reese, Rincon, and Rieu- 
tord (2004, RRR) have found the eigenmodes of a shell of a 
conductive incompressible fluid with a dipole magnetic field, 
[see also Rincon and Rieutord 2003, and Levin and D'Angelo 
(2004)]. However, we are not just limited by purely computa- 
tional difficulties. In a magnetar core the magnetic pressure 
is sub-dominant, and there is virtually no conversion of the 
Alfven waves into magnetosonic waves. The core responds to 
the wave as an incompressible fluid, and the wave propagates 
along the field lines. Since different field lines have different 
lengths in the core, and different average Alfven velocities 
along them, one might suspect that the core supports a con- 
tinuum of the Alfven modes. This is hard to prove generally, 
but below we, in the spirit of Appendix B in RRR, present 
a toy-model example where the Alfven-mode spectrum can 
be computed explicitly. 

Consider a perfectly conducting incompressible fluid 
sandwiched in a box, with top and bottom plates also being 
perfect conductors; see Figure 1. The magnetic field B(y) 
is vertically directed and is a function of y only, threading 
both top and bottom plates. Gravity is zero, and magnetic 
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pressure gradient is trivially compensated by the fluid pres- 
sure gradient. Consider now a ^-independent Lagrangian dis- 
placement y) in the z direction. The equation of motion 
is 



dt 2 



c 2 (y) 



2 C 



9a; 2 



■7 



dt' 



(4) 



where c(y) — \jT{y)j p(y) and 7 is a small damping con- 
stant added for generality; here T = B 2 /(4n) is the mag- 
netic tension (in a superconducting fluid T = £?B C rit/(47r)) 
and we have allowed y-dependence of the fluid density p. 
Clearly, Eq. (4) separates in x and y, with eigenfunctions 



Cny (x,y) = sm(nnx/l x )S(y - y ) exp(iuj nya t) , 
with eigenfrequencies 



wny = ±V( n */ l *) 2T (y°)/p(.y°) - 774 + »7/2- 



(5) 



(6) 



Here l x is the box height, n is the integer number, and yo is 
the y-coordinate where the eigenfunction is localized. The 
spectrum is continuous! 

It is notable that for a spherical shell numerical results 
of RRR indicate that the spectrum of toroidal modes is also 
continuous [and the eigenmodes are singular, as in Eq. (5) 
above]. The continuity of part of the spectrum is a topologi- 
cal property, which should remain when external parameters 
(the shape of the box, density profile, etc) change continu- 
ously. Since the basic MHD physics in our example, in the 
fluid shell of RRR, and in the magnetar core are similar, 
we believe that the magnetar core supports a continuum of 
MHD modes. 

We started this section by asking whether a global mag- 
netar MHD-elastic mode can be responsible for an observed 
QPO. We can now see that such a mode of the whole mag- 
netar will be coupled to a continuum of the MHD modes in 
the core. Experience from Quantum Mechanics tells us that 
a mode coupled to a continuum of other modes decays if 
its frequency lies within the continuum (cf. Fermi's Golden 
Rule). This is likely the situation here: the QPO frequencies 
are tens of Hertz, above the base AlfVen-mode frequency 
(but obviously below that of infinite number of higher-order 
Alfven modes). Let us illustrate explicitly the decay of a 
global mode by returning to the simple example in Figure 
1. We now model the presence of a crust by assuming that 
the top plate is allowed to move in z-direction but is also a 
harmonic oscillator (with restoring force provided by some 
external spring) . The equation of motion of the top plate is 
given by 



(7) 



where Z is the top plate's displacement, uio is the proper 
frequency of the oscillator, F is the force due to magnetic 
stresses from the fluid below, and M is the plate's mass. 
When the top plate participates in a global mode of fre- 
quency w, then Z = Zoexp(iiot). It is straightforward, by 
using Equation (4), to solve for the fluid motion: 



S m[k(y)x] 
^ X ' V) S m[k(y)l x ] Z ' 

where 

k{y) = \/u! 2 - iu>-y/c{y) 
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(10) 



The back-reaction force F can now be computed: 
F = -/T(,)(g)^ 

= J T(y)k(y)cot[k(y)l x ]dy 

where l z is the z-dimension of the top plate. Let 7 = 0. 
Then if to real and is in resonance with at least one of the 
continuum modes then cot kl x will have a pole singularity at 
some y = yo in the range of integration. The singularity is 
lifted if 7 is non-zero; in the limit of 7 -C w the integration 
will produce an imaginary component of the amplitude of 
F: 



Q[F exp(-iwt)] = -IzCo 



T(y )k(y )n 



(11) 



Comparing the above equation with the Eq. (7), we see 
that our assumption that to was real was not self-consistent. 
Within an order-of-magnitude the imaginary part of ui is 



l x Mtu 

Substituting the parameters relevant for a magnetar (B 
1Q 15 i z — 1 1 — m6™ in3l 

we get 

~ 10s 



(12) 



g, u> ~ 300rad/s), 
(13) 



Even allowing for several corrections of order 1, we see that 
a global mode is expected to decay within a second. This 
is much shorter than the timescale on which the QPO is 
observed, ~ 100s. 

One can ask the following question: is it possible that 
a global mode is not coupled to the MHD continuum in the 
core? After all, along with the MHD continuum the core 
can probably support discrete MHD modes as well (such 
modes are seen in RRR's numerical work). However, while 
some global mode may be decoupled from the continuum, 
the mode of our interest must involve the motion of the 
crust, in order to produce the observed QPO. Unless mag- 
netic field in the core has high degree of symmetry (which 
seems unlikely, especially if one believes that it was gener- 
ated by turbulent dynamo), this crustal motion is bound to 
couple the QPO-generating global mode to the continuum 
of the core's MHD modes via the field-lines threading the 
crust. Thus the decay of a global mode seems generic. 



4 MAGNETAR: A TUNING FORK OR WET 
SPAGHETTI? 

Observations argue for the former: QPO frequencies are sta- 
ble over tens of seconds, something strongly indicative of a 
mechanical mode. However, we have shown that the theory 
argues for the latter. Strong MHD crust-core coupling de- 
stroys the stability of a purely crustal mode, and the pres- 
ence of a continuum of MHD modes in the core makes a 
global crust-core mode decay within a second. So what is 
the origin of the magnetar QPOs? It can still be mechani- 
cal, but only if the magnetic field geometry is very different 
from what has been assumed — for example, if the magnetic 
field is largely confined to the crust or has an extremely in- 
coherent poloidal component. Alternatively, the QPOs may 
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have a magnetospheric origin; this was mentioned as a pos- 
sibility in Barat et al. (1983), but has remained unexplored. 
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